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AN  ANALYSIS  OF  THE  FREE  VIBRATION  OF  A SHALLOW 
SPHERICAL  MEMBRANE  SHELL 


INTRODUCTION 

The  ferroelectric  element  in  some  transducers  is  a shell  that  has  the 
form  of  a sector  of  a sphere.  In  order  to  design  practical  underwater 
acoustic  transducers  that  use  spherical-shell-sector  elements,  it  is 
useful  to  know  the  lowest  resonant  frequencies  of  such  elements  and  to 
establish  the  modal  displacement  functions  that  correspond  to  these 
resonant  frequencies.  This  report  describes  the  calculation  of  the 
lowest  resonant  frequencies  and  the  corresponding  modal  displacement 
functions  of  a shallow  spherical  membrane  shell  that  is  composed  of 
a homogenous  material.  The  analysis  presented  here  does  not  provide 
a complete  description  of  the  vibration  of  a piezoceramic  transducer 
element  of  this  form,  in  that,  it  does  not  explicitly  account  for  the 
effect  of  the  electrical  boundary  conditions  upon  the  motion  of  the 
shell.  Nevertheless,  this  analysis  is  necessary  before  design 
procedures  for  transducers,  which  incorporate  piezoelectric  shells  that 
are  shallow  sectors  of  spheres,  can  be  established.  A comprehensive 
theory  that  treats  a spherical  shell  composed  of  a ferroelectric 
material,  and  which  incorporates  the  complete  set  of  electromechanical 
equations,  would  be  the  next  step  in  the  analytic  treatment  initiated 
here. 


Not*:  Manuscript  submitted  July  IT.  1978. 


The  subsequent  material  in  this  report  will  be  presented  as  follows. 
First,  a description  of  the  spherical  shell  will  be  given,  and  those 
formulas  from  differential  geometry  that  will  be  needed  in  the  later 
analysis  will  be  written  down.  Also, the  assumptions  and  limits  of 
shallow-shell  theory  will  be  noted.  Next,  the  equations  of  motion  of 
a shallow  spherical  membrane  shell  will  be  written  down  in  several 
different  forms.  The  equations  will  Le  solved  for  the  normal  and 
tangential  displacement  functions.  Permissible  boundary  conditions 
will  be  discussed.  The  corresponding  secular  equations,  whose  roots 
yield  the  eigenf requencies  of  the  shell's  free  vibration,  will  be  given. 
In  addition  to  these  secular  equations,  which  result  from  the  solution 
to  the  equations  of  motion  that  describe  the  shell,  an  approximate 
secular  equation  will  be  derived  by  use  of  the  Ray leigh-Ritz  method. 

The  analysis  in  this  report  is,  for  the  most  part,  a codification  of 
a number  of  the  results  found  in  the  excellent  book  on  shell  theory  by 
Kraus  [l]  , and  it  draws  freely  and  extensively  upon  the  material  in 
this  text.  It  is  desirable,  however,  to  collect  Kraus'  results  for 
shallow  spherical  membrane  shells  and  to  present  them  in  a form  that  is 
useful  in  transducer-design  work.  Moreover,  several  errors  were  found 
in  rederiving  Kraus'  equations.  These  have  been  corrected  in  the 
equations  reported  here.  Those  interested  in  further  details  of  the 
analysis  presented  in  this  report  should  refer  to  the  comprehensive 
treatment  in  Kraus'  text.  The  notation  used  in  this  report  is  the  same 
as  that  used  by  Kraus,  whenever  possible. 

GEOMETRY  OF  A THIN  SHALLOW  SPHERICAL  SHELL 

Since  a thin  shell  can  be  considered  as  the  materialization  of  a 
curved  surface,  results  from  differential  geometry,  which  is  needed  to 
describe  such  surfaces  in  space,  are  incorporated  into  the  theory  of 
thin  shells.  The  motion  of  a point  in  a thin  shell  is,  in  fact, 
determined  by  the  motion  of  the  curved  reference  surface  associated 
with  the  shell.  Consider  the  geometry  of  a rotationally  symmetric  thin 
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shell  of  general  form  that  Is  shown  in  Fig.  1.  The  reference  surface  S 
of  this  shell  may  be  described  in  the  rectangular  coordinate  system 


SHELL 
REFERENCE 
SURFACE  S 


MERIDIONAL  CURVE 
(GENERATING  CURVE) 


LATITUDE 

CIRCLE 


Fig.  1 - Geometry  of  a rotationally  symmetric  thin  shell 


(Xj.x^.x^)  as  a function  of  two  independent  curvilinear  coordinates  of 
the  surface,  <*.  and  o , by  the  expressions 


The  surface  S is  symmetric  with  respect  to  the  x.  axis.  The  first 

J ■> 

fundamental  form  of  the  surface  S is  a relation  that  describes  (dr)“, 
the  square  of  the  length  of  the  infinitesimal  vector  dr  that  measures 


I 


Che  change  in  Che  vecCor  r as  one  moves  from  poinC  P Co  an  infiniCesi- 
mally  near  neighboring  poinC  on  Che  surface.  If  and  are  orchog- 
onal  curvilinear  coordinaces,  Chen  Che  firsC  fundamencal  form  of  S is 


(ds)2  - A2(d0l)2  + A2(dc2)2 


(2) 


in  which  A^  and  A2  are  fundamencal  consCanCs  characCerizing  S.  The 
firsc  fundamencal  form  can  also  be  wriccen 


(ds)2  - R2(d<fr)2  + R2(d9) 2 


(3) 


in  Cerms  of  Che  independenc  coordinaces  $ and  9,  where  41  is  Che  meridi- 
onal angle,  ChaC  is,  Che  angle  beCween  Che  shell's  axis  of  rocacional 
symmecry  and  Che  normal  Co  S at  P and  where  9 is  Che  azlmuChal  angle  of 
Che  meridional  plane  concaining  P.  The  quanCiCy  R^  in  Eq.  (3)  is  Che 
principal  radius  of  curvaCure  of  Che  meridional  curve,  which  [2]  lies 
along,  buC  has  a direcCion  opposiCe  Co,  Che  normal  Co  Che  surface  S ac 
P.  The  quanCiCy  RQ  is  Che  radius  of  Che  laClCude  circle  ChaC  passes 
Chrough  P.  Comparing  Eqs.  (2)  and  (3),  one  sees  ChaC  a^  and  may  be 
associaced  wich  41  and  9,  respecCively , and  ChaC 


and 


A„  - R„ 


(4a) 


(4b) 


Three  differencial  equaClons,  Che  Gauss-Codazzl  condlCions,  relace  Che 
fundamencal  consCanCs  A^  and  A^  Co  Che  principal  radii  of  curvacure  of 
a surface  at  a poinC.  When  applied  in  Che  (4,9)  coordinace  syscem,  Che 
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Gauss-Codazzi  conditions  require  that 


dR0 

jf  " Rlcos^  (5> 

The  foregoing  concepts,  which  apply  to  a general,  rotationally 
symmetric  shell,  can  now  be  specialized  to  a shallow  spherical  shell. 
A shell  that  is  the  surface  of  a segment  of  a sphere  is  depicted  in 
Fig.  2.  The  principal  radius  of  curvature  R^  is  constant  for  such  a 
spherical  shell.  This  radius  of  curvature  has  been  denoted  as  R, 
rather  than  R^,  in  Fig.  2.  It  is  easily  seen  that  the  origin  of  the 
(x^.x^.x^)  system  of  Fig.  1 may  be  placed  so  that  R terminates  at  the 


Fig.  2 - Geometry  of  a spherical  shell 


center  of  that  sphere  from  which  the  shell  is  derived.  Also,  for  this 
spherical  shell,  the  meridional  radius  has  been  denoted  by  r instead  of 
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by  Rq.  Thus,  in  Che  revised  notation,  Eqs.  (4)  become 

Aj  - Rj  - R (6a) 

and 

A 2 ■ Rq  “ r.  (6b) 

Figure  2 shows  a to  be  the  radius  of  the  spherical  shell  at  its  rim, 
to  be  the  meridional  angle  of  points  on  the  rim,  and  H to  be  the 
height  of  the  shell  segment.  It  is  also  readily  shown  that,  for  a 
general  spherical  shell, 

aJ  - 2RH  - H2.  (7) 


The  notion  of  shallowness  is  introduced  into  thin-shell  theory  in 
order  to  simplify  the  equations  that  govern  the  motion  of  a shell. 
Reissner's  [3]  criterion  for  considering  a shell  shallow  is  that 


which  implies  that  a spherical  shell  may  be  considered  shallow  if 

i 30°.  (10b) 

If  $ is  small,  one  has,  from  Eq.  (5),  the  result  that 

dr  ss  Rd;  <ll) 

for  a spherical  shell.  Equation  (11)  is  basic  in  the  theory  of 
shallow  shells,  even  though  $ is  not  small  for  all  shells  that 
satisfy  Reissner's  criterion  of  shallowness;  i.e.,  all  those  for 
which  Eq.  (10b)  holds. 

FREE  VIBRATION  OF  A SHALLOW  SPHERICAL  MEMBRANE  SHELL 

Thin-shell  theory  is  usually  considered  valid  if  the  thickness 
h of  the  shell  in  question  is  small  in  comparison  to  one  of  its 
radii  of  curvature.  A spherical  shell  is,  thus,  thin  if 

h/R  < £,  (12) 

where  e in  engineering  studies  [•*]  is  usually  taken  as  1/10  or  1/20. 
Under  conditions  where  a thin  shell  is  so  loaded  that  all  bending 
moments  are  zero  or  are  negligibly  small,  the  shell  is  said  to  be 
in  a membrane  state  of  stress.  The  set  of  equations  that  governs 
the  motion  of  a thin  shell  is  considerably  simplified  if  the  shell 
is  in  a state  of  membrane  stress,  in  which  case  it  is  referred  to 
as  a membrane  shell.  The  membrane-shell  assumptions  will  be  considered 
to  hold  for  the  spherical  shells  analyzed  subsequently. 
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The  equations  of  motion  of  a shallow  spherical  shell  of  thickness 
h,  which  is  not  loaded  by  external  forces  and  which  is  executing  motion 
that  is  independent  of  the  circumferential  position,  i.e.,  independent 
of  the  coordinate  0,  are 


3(rN  ) 

— Ne 


» prh 


2 

u 

9* 

3t ' 


(13a) 


3(rQ  ) .2 

_1  . I(Nr  + Nd)  - orh  i-|. 


(13b) 


and 


3(rM  ) 

“IT  ‘ M0  - rQr  * °* 


(13c) 


in  which  p denotes  the  density  of  the  shell  material,  u and  w, 

respectively,  denote  the  tangential  and  the  normal  shell  displacements, 

N and  N,  are  the  stress  resultants,  M and  M,  are  the  bending  moments, 
r6  r 9 6 

is  the  shear-force  resultant,  t denotes  the  time,  and  r and  R are 
as  shown  in  Fig.  2.  For  a membrane  shell, 

M - M « 0,  (1‘a) 

r 9 

and,  consequently,  by  Eq.  (13c) 


Qr  - 0.  (14b) 

As  a consequence  of  Eqs.  (14),  the  three  equations  of  motion  reduce  to 
two  equations,  Thus,  if  the  forces  and  displacements  are  assumed  to 
have  a harmonic  time  dependence  described  by  the  factor  exp(jwt) 

(which  is  henceforth  dropped),  Eqs.  (13)  become,  for  a shallow  spherical 
membrane  shell  in  free  harmonic  vibration,  the  two  coupled  differential 
equations 
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1 


(15a) 


Figure  3 shows  a differential  element  of  the  spherical  shell  in 
order  to  illustrate  the  action  of  the  forces  N and  N,  and  the 


SHELL 

ELEMENT 


Fig.  3 - Differential  element  of  a shallow  spherical  membrane  shell 


directions  of  the  displacement  u and  w.  The  shell  element  shown  in 
Fig.  3 is  defined  by  two  infinitesimally  separated  meridional  curves 
and  by  two  Infinitesimally  separated  latitude  circles.  The  stress 
resultant  N acts  normal  to  the  face  abed  of  the  element,  and  the  stress 
resultant  N acts  normal  to  the  element  face  edef.  The  direction  of 
the  displacement  w is  normal  to  the  reference  surface  of  the  shell, 
that  is,  normal  to  element  face  befg,  and  the  direction  of  the 
displacement  u is  in  the  meridional  direction.  Owing  to  the  shell's 
symmetry,  there  can,  of  course,  be  no  displacement  in  the  d direction. 
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If  the  results  expressed  by  Eqs.  (16)  through  (18)  are  substituted 
into  Eqs.  (15),  one  obtains,  for  the  two  equations  of  motion,  the 
expressions 


L(ru)  + w'  [r ( 1 + v)/r]  + (phui  ru/K)  ■ 0 


and 


(ru)'  + (2rw/R)  - |[phRw2w]/  [K(l  + v)  ] \ • 0. 

In  Eq.  (20),  the  prime  denotes  differentiation  with  respect  to  r 


(19) 


(20) 


(•••)'  - (...), 


and  the  operator  L in  Eq.  (19)  is  defined 


L( . . . ) = (...)"  - ~ (...)' 


(21a) 


d “ ^ d L 
r dr2  dr  “ r 


(...). 


(21b) 


The  differential  equations  for  the  normal  and  tangential  displacements 
of  the  shell,  w(r)  and  u(r),  are  coupled  when  expressed  in  the  form  of 
Eqs.  (19)  and  (20).  Uncoupled  differential  equations  for  w and  u may, 
however,  be  found  in  the  following  way.  Suppose  one  solves  Eq.  (20) 
for  w,  and  obtains  the  result 


w - [u1  + (u/r)  ]/|3. 


(22) 
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where 


8 - (phu)2R)/[K(l  + v)]  - (2/R)  • (23) 

Then,  if  Eq.  (22)  Is  differentiated  with  respect  to  r,  the  expression 

w*  - L(ru)/(8r)  (24) 

Is  produced.  Substituting  the  value  of  w' , given  by  Eq.  (24),  into 
Eq.  (19)  allows  one  to  obtain  the  uncoupled  differential  equation  for 
the  tangential  displacement  of  the  shell,  viz 

L(ru)  + (p/a)“(ru)  - 0,  (25) 


where 

u2  - a2tf2(l  + v)[2  - tf2(l  - v)  ]/[R2(l  - 02)], 

in  which 

2 2 ^ 

« - pw  R /E. 

Now  Eq.  (19)  can  be  solved  for  u,  with  the  result  that 


(2b) 


(27) 


u - -[K/(orhio")  ] { [r(l  + v)/R]  w'  + L(ru)}.  (28) 

Upon  substituting  the  value  of  L(ru)  that  is  obtained  from  Eq . (24)  into 
Eq.  (28),  one  obtains  the  result  that 


- [R(  1 - a2)  ]/ [iJ2 ( 1 + v)]w\ 


U 


(29) 


t 


which  yields,  upon  differentiation,  the  equation 

u'  - [R(l  - fl2)]/[rt2(l  + v)]  w".  (30) 

The  expressions  for  u and  u' , given  by  Eqs.  (29)  and  (30),  when 
substituted  into  Eq.  (22),  results  in  the  uncoupled  differential 
equation  for  the  normal  displacement  of  the  shell,  viz 

w"  + w'/r  + (p/a)“w  » 0.  (31) 

When  one  differentiates  Eq.  (31)  with  respect  to  r and  multiplies  both 
sides  by  r,  he  obtains  the  equation 

L(rw')  + (u/a)“(rw')  “ 0.  (32) 

Since  Eqs.  (25)  and  (32)  have  the  same  form,  their  respective  solutions 
u(r)  and  w'(r)  must  be  proportional.  This  can  easily  be  shown  by 
substituting  the  result 

L(ru)  - -(u/a)“ru,  (33) 

obtained  from  Eq.  (25),  into  Eq.  (19).  When  this  is  done,  one  obtains 
the  result  that 


u/V  - R(1  + v) / [( pR/ a) 2 - (1  - v2)fl2].  (34) 


One  notes  that  if  the  change  of  variable  s - (p/a)r  is  made  in 
Gq.  (31),  the  result  is  Just  Bessel's  equation  of  order  zero.  The 
general  solution  to  Eq.  (31)  therefore  is 

v - -A(a/p) Jq( pr/a) , (35) 

where  Jrt(...)  is  the  Bessel  function  of  the  first  kind  of  order  zero 
and  where  A is  an  arbitrary  constant.  The  zero-order  Bessel  function 
of  the  second  kind  Y^(...)  does  not  enter  the  general  solution  to 
Eq.  (31),  since  w(r)  cannot  be  singular  at  r - 0.  From  Eqs.  (34)  and 
(35),  one  also  obtains  the  result 

u - {R(l  + v)  / [( pR/ a)  2 - (1  - v2)il2]}  AJ^pr/a),  (3b) 

in  which  J^(...)  is  the  Bessel  function  of  the  first  kind  of  order  one. 

Either  of  two  possible  boundary  conditions  may  be  imposed  at  the 
edge  r - a of  a shallow  spherical  membrane  shell.  At  r ■ a,  one  can 
have  either 

(37) 

(38) 

The  boundary  condition  that  is  expressed  by  Eq.  (37)  corresponds  to 
a tangentially  clamped  shell  rim,  while  that  expressed  by  Eq.  (38) 
corresponds  to  a free  edge.  Note  that  neither  of  the  two  permissible 
boundary  conditions  is  imposed  upon  the  normal  displacement  w. 


or  else 


u(a)  ■ 0 


Nr(a) 


- 0. 


Equation  (36)  requires  that  y must  be  a root  of  the  secular  equation 


JL(h)  - 0 


for  the  tangentially-clamped  boundary  condition  that  is  expressed  by 
Eq.  (37).  For  the  boundary  condition  expressed  by  Eq.  (38),  which 
characterizes  a shell  with  a free  edge,  one  finds  that  y must  be  a root 
of  the  secular  equation 


JQ(u)  - |(uR2/a2)/[(yR/a)2  - fl2(l  - v2)]} 

X [ ( v - l)Jj(y)  + yJg( y)  ] - 0.  (AO) 

Equation  (40)  is  obtained  as  follows.  First,  one  calculates  the  various 
terms  appearing  in  the  stress  resultant  N , which  is  given  by  Eqs.  (16a) 
and  (18),  using  values  for  u,  u'  and  w obtained  from  Eqs.  (35)  and  (36). 
One  then  sets  the  result  obtained  equal  to  zero. 


The  roots  of  Eq.  (40)  determine  the  resonant  frequencies  of  the 

spherical  shell  with  a free  edge.  Before  Eq.  (40)  can  be  solved, 

2 

however,  it  is  necessary  to  express  the  quantity  0 therein  in  terms 
of  u.  To  do  this, one  rewrites  Eq.  (26)  in  the  form 


04  - {[2(1  + v)  + (yR/a)2]/ (1  - v2) } fl2  + (yR/a)2/(l  - v2)  - 0 (41) 


2 2 

and  solves  this  quadratic  equation  in  0 for  its  two  roots  fl.(y)  and 

2 2 1 
rt„(y).  To  each  of  these  two  values  of  0 there  is  a set  of  roots  u of 
2 2 

Eq.  (40).  When  inserted  into  Eq.  (40),  the  root  (y)  generates  a 
sequence  of  eigenvalues  u^,  with  {m  * 1,  2,  •••}  , when  that  equation  is 

solved.  Another  sequence  of  eigenvalues  y.  , again  with  (m  ■ 1,  2,  ...l, 

2m  ‘ 


2 

is  obtained  when  this  procedure  is  repeated  with  ft7(y).  There  are, 
thus,  two  families  of  solutions  y^,  with  k ■ 1 or  k * 2,  to  Eq.  (40). 

To  obtain  the  resonant  frequencies  of  the  shallow  spherical  membrane 
shell,  each  root  y^  is  inserted  into  Eq.  (41).  Each  of  the  two  values 
of  y^  yields  in  turn  two  values  of  ft^ , which  can  be  designated 
2 

and  There  thus  appear  to  be  four  fiimilies  of  resonant  frequencies 

of  the  shallow  spherical  membrane  shell  with  a free  edge.  Recalling 
Eq.  (27),  one  notes  that  the  resonant  angular  frequencies  are 


ikm 


«W»  (E/o) 


(42) 


with  i - 1 or  i ■ 2 and  k - 1 or  k ■ 2.  Of  course,  only  those  roots 

y^  which  yield  real  values  of  the  resonant  frequencies  can  be 

considered  to  be  of  physical  significance.  (One  notes  in  passing  that 

there  are  only  two  families  of  resonances  for  a spherical  shell  with  a 

2 

tangentially  clamped  edge.  Since  Eq.  (39)  does  not  involve  0 , it 

yields  directly  a sequence  of  eigenvalues  y . Each  of  these  y when 

m , m 

Inserted  into  Eq.  (41)  gives  rise  to  two  values  of  ft  , which  can  be 
2 2 

denoted  ftm1  and  ft^.,.  One  resonant  frequency  co^  of  the  shell 
corresponds  to  each  such  ft^-) 


The  modal  displacement  functions  corresponding  to  the  roots  of 

Eq.  (40)  can  now  be  expressed.  First,  note  from  Eq . (35)  that  the 

expression  for  the  normal  displacement  w does  not  depend  explicitly 

upon  ft.  Thus,  there  are  just  two  families  of  modal  normal  displacement 

functions  corresponding  to  the  two  sequences  of  eigenvalues  y,  and 

lm 

y_  . These  modal  normal  displacement  functions  are 


w 


km 


-(a/ukm)J0(\mr/a)‘ 


(43) 
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with  k ■ 1 or  k ■ 2.  On  the  other  hand,  since  by  Eq.  (36)  the 
tangential  displacement  depends  upon  ft,  there  appear  to  be  four 
families  of  modal  tangential  displacement  functions,  namely 

“ik.  ■ lR<‘  * - u - <**> 

with  i ■ l or  i ■ 2 and  k ■ 1 or  k ■ 2.  The  modal  functions  of  the 
gravest  modes  of  the  shallow  spherical  membrane  shell  are  those 
corresponding  to  setting  m equal  to  one  in  Eqs.  (43)  and  (44). 

I ,| 

MODAL  VIBRATIONS  OF  A SHALLOW  SPHERICAL  MEMBRANE  SHELL  WITH  A FREE 
EDGE,  CALCULATED  USING  THE  RAYLEIGH-RITZ  METHOD 

In  order  to  calculate  the  resonant  frequencies  of  free  vibration 

of  a shallow  spherical  membrane  shell  that  has  a free  edge,  one  is 

required  to  find  the  roots  of  Eq.  (40) . Moreover,  as  was  explained 

previously,  before  Eq.  (40)  can  be  solved,  it  is  necessary  first  to 

2 

solve  Eq.  (41)  for  jl  in  terms  of  u and  to  insert  each  of  the  two 
2 

values  of  £1  thus  obtained  into  Eq.  (40)  before  proceeding  to  find  the 
roots  of  this  latter  equation.  While  such  a procedure  for  calculating 
the  resonant  frequencies  of  a spherical  shell  with  a free  edge  is 
clearly  possible,  it  is  sufficiently  complicated  to  make  simpler 
approximate  methods  of  calculating  shell's  resonant  frequencies 
attractive. 

[ i 

One  such  approximate  method  of  calculating  the  resonant  frequencies 
of  a shell's  vibration  is  the  Rayleigh-Ritz  method.  The  formulation  of 
the  Rayleigh-Ritz  analysis  that  is  appropriate  in  thin-shell  theory 


ft 


i 
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is  based  upon  Che  variacional  equation 


{[  ^^(u.v.w)  + phuTu]6u 


+ [ 2(u,v,w)  + phoTv]6v 


+ [ + phu^wJdw^A.,  dc^do^  - 0. 


(45) 


in  which  the 


) are  the  differential  operators  that  appear  in 


the  equations  of  motion  for  free  vibration  of  a thin  shell,  once  the 
shearing  forces,  Che  stress  resultants,  and  the  bending  moments  have 
been  eliminated.  In  terms  of  these  operators,  the  equations  of  motion 
appear  in  the  form 


J^\(u,v,w)  - ph  , 
it 


(46a) 


(u,v,w)  - ph  , 


(46b) 


3t 


and 


(u,V,w)  - ph 

3t 


(46c) 


For  symmetric  free  vibration  of  the  kind  being  considered,  there  is  no 
tangential  displacement  of  the  shell  in  the  direction  of  9.  Therefore, 
for  symmetric  motion,  the  displacement  v does  not  appear  in  Eq.  (45), 
and  one  has  a variational  equation  of  the  form 
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Li L 


U ^(.u.w)  + phu)‘u]6u 


ai  a: 


+ [ 5^,(u,w)  + phui'wjiwlAjA^  da^da.,  - 0. 


Likewise  for  symmetric  motion,  there  is  no  appearance  of  v in  Eqs.  i4b) 


For  the  shallow  spherical  shell,  one  has,  from  the  previous 
discussion  of  shell  geometry  and  from  Eqs.  (b)  and  (11),  the  result 


A^A.,  do^da,,  • r drdi>. 


Since  the  shell's  motion  is  independent  of  t),  the  integral  over 

a.,  ■ S can  be  carried  out  directly  after  one  substitutes  the  result 

of  Eq.  (48)  into  Eq . (47).  The  resulting  expression  is  a variational 


integral  over  r: 


[1%  (u,w)  + phu)"u]5u 


♦ tiR  (u,w)  + ohui-wl^fr  dr  - 0 . 

* 

Note  that,  in  the  case  of  harmonic  free  vibration  of  a shallow 
spherical  membrane  shell,  Eqs.  (19)  and  (20)  can  be  put  in  the  form 
of  the  general  equations  of  shell  motion  that  are  expressed 
symbolically  by  Eq.  (46).  When  this  is  done,  one  can  obtain  from 
Eq.  (49),  for  the  shell  in  question,  the  expression 


L 
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- C{(1  + v) /R} { (ru )'  + (2rw/R)} 

- ( iJ/ R)  *~ ( 1 - v“)rw]6wj>  dr  - 0.  (50) 


In  obtaining  Eq.  (50),  the  results  expressed  by  Eqs.  (17)  and  (27) 
have  been  used.  Equation  (50)  is  the  form  of  the  Ravleigh-Ritz 
variational  equation  that  is  needed  for  finding  the  resonant 
frequencies  of  a shallow  spherical  membrane  shell. 

The  next  step  in  the  Rayleigh-Ritz  method  is  to  find  suitable 
approximate  displacement  functions,  that  is,  suitable  approximate 
expressions  for  u and  w.  In  order  to  use  Eq.  (50),  however,  these 
approximate  displacement  functions  must  be  such  that  the  boundary 
conditions  appropriate  to  the  problem  in  question  are  satisfied. 

For  the  vibrating  spherical  shell  with  a free  edge,  this  boundary 
condition  is  given  by  Eq.  (38),  which,  when  written  in  terms  of 
the  displacement  functions  u and  w,  has  the  form 

u'  + [(l  + v)/R]w  + (v/r)u  - 0.  (51) 

Equation  (51)  holds  at  r • a. 
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Approximate  forms  of  the  displacement  functions  are  found  in  the 
following  way.  First,  note  that  the  exact  forms  of  the  displacement 
functions  are  available.  In  particular,  one  has  from  Eq.  (3o)  that 


u - BJ1 (yr/a) , 


where  B is  a constant.  Consider,  on  the  one  hand,  the  form  of  the 
series  expansion  of  a first-order  Bessel  function: 


w / \ x , x ^ x . 

Jl(x)  " 2 ( 1 " 8 + 19 2 - } • 


(53a' 


Now  suppose,  on  the  other  hand,  that  one  makes  the  approximation 


2 2 .4  4 

J^(x)  a sin(bx)  - bxt  I + 


(53b) 


in  which  b is  a constant  chosen  so  as  to  make  the  first  zero  of  J^(x) 
coincide  with  the  first  axis-crossing  of  the  sine  function.  The 
approximation  expressed  by  Eq.  (53b)  can  be  expected  to  be  reasonably 
good  for  values  of  x somewhat  beyond  the  first  zero  of  J^(x).  Upon 
examining  Eqs.  (52)  and  (53),  it  is  not  too  difficult  to  see  that  a 
good  approximation  of  the  displacement  function  u would  be 


u - A(. r/a)  [l  - Y(r/a)-],  (54) 

in  which  A is  an  arbitrary  constant  and  in  which  > is  a constant  chosen 
in  such  a way  that  the  boundary  condition  in  question  is  satisfied. 

For  example,  in  the  case  of  the  tangentially  clamped  shell, for  which 
the  boundary  condition  is  expressed  by  Eq.  (37),  one  would  take  y - 1. 
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The  approximate  displacement  function  expressed  by  Eq . (54)  can  be 
expected  to  be  a good  approximation  of  the  tangential  modal 
functions  corresponding  to  the  lowest  resonant  frequencies, 

Since  the  boundary  condition  for  a spherical  shell  with  a 
free  edge  is  expressed  by  Eq.  (51),  an  approximation  for  the 
displacement  function  w,  corresponding  to  the  approximation  for 
u that  is  expressed  by  Eq.  (54),  is  also  needed.  This  may  be  readily 
found  by  inserting  the  value  of  u expressed  by  Eq . (54)  into  Eq.  (22). 
The  result  is  that 

w - A[2/(afl)][l  - 2y(r/a)2],  (55) 

in  which,  after  combining  Eqs.  (17),  (23),  and  (27),  one  expresses  g 
by  the  equation 

0 - [(1  - v)fl2  - 2]/R.  (56) 

Note  that,  since  Eq.  (22)  is  just  a rearrangement  of  Eq.  (20),  the 
approximate  displacement  functions,  given  by  Eqs.  (54)  and  (55),  will 
satisfy  Eq.  (20). 

The  approximate  displacement  functions  u and  w,  which  are  given  by 
Eqs.  (54)  and  (55),  are  now  substituted  into  Eq.  (51),  which  expresses 
the  boundary  condition  for  the  shell  with  a free  edge.  The  resulting 
equation,  when  solved  for  y,  yields  the  expression 

Y - [(1  + v)Q2]/[(3  + v)02  - 2].  (57) 
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Thus,  if  y is  given  the  value  expressed  by  Eq.  (57),  the  approximate 

displacement  functions  will  satisfy  the  proper  boundary  condition. 

2 

Note  from  the  presence  of  0 in  Eq.  (57)  that  both  of  the  approximate 
displacement  functions  u and  w must  depend  upon  the  frequency  of 
vibration  of  the  shell,  if  the  rim  of  the  vibrating  shell  is  to  be 
stress  free. 


Equation  (50)  can  now  be  used  to  derive  an  approximate  secular 
equation  for  the  resonant  frequencies  of  a shallow  spherical  membrane 
shell.  First,  note  that  from  Eq.  (54)  one  has,  for  the  variation  in  u, 
the  expression 


6u  » (r/a)[l  - y(r/a)^]6A.  (58) 

Next,  the  displacement  functions  given  by  Eqs.  (54)  and  (55)  and  the 
result  expressed  by  Eq.  (58)  are  substituted  into  Eq.  (50) . Because 
the  displacement  functions  u and  w satisfy  Eq.  (20),  the  bracketed 
term  in  the  integral  that  multiplies  6w  will  vanish  identically. 
Therefore,  one  need  not  deal  further  with  that  term.  The  substitution 
in  question,  which  is  somewhat  simplified  if  one  recalls  Eq.  (24), 
yields  the  equation 


■’I'-’*) 


(59) 


When  the  integration  in  Eq.  (59)  is  performed,  the  equation 
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A6A{y(1  - ir)(3  - 2y)  / [( 1 - v)iT  - 2] 

+ + v)(3y“  - 8y  + 6)/(4R)~]}  - 0 

Is  obtained.  Since  the  variation  <$A  is  arbitrary  and  since  the 
constant  A must  be  non-zero  for  non-trivial  displacement  functions, 
Eq.  (60)  is  satisfied  only  if  the  term  within  braces  vanishes.  If 
the  value  of  y given  by  Eq.  (57)  is  used,  Eq.  (60)  with  the  term 
in  braces  set  equal  to  zero  leads  to  the  following  secular  equation: 


c.o6  + + c,o2  + C - 0, 

l Z J 


in  which  the  constants  are 


- (l  - v) (33  + lOv  + v~) , 


(62a) 


C2  --[2(61  - 14 v - 3v2)  + (4R/a)2(7  + v)], 


(62b) 


C3  - 8(17  - v)  + (4R/a)  *■(  13  + v) , 


(62c) 


C4  - -48[2(R/a)  + 1 ]. 


(62d) 


2 2 

Equation  (61)  is  cubic  in  £ . Suppose  is  the  smallest  real 
root  of  Eq.  (61).  The  resonant  angular  frequency  of  the  gravest  mode 
of  a shallow  spherical  membrane  shell  with  a free  edge  is,  therefore. 


ui0  - (00/R)(E/p)\ 
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It  Is  known  that  the  Rayleigh-Ritz  method  always  yields  a resonant 
frequency  that  Is  somewhat  higher  than  the  exact  result  (i.e.,  the 
resonant  frequency  determined  from  Eq.  (AO).  Therefore,  the 
approximate  resonant  frequency  of  a shallow  spherical  rnembpame  shell 
given  bv  Eq.  (63)  will  be  somewhat  greater  than  tJia*'"found  using  the 

exact  theory  that  was  outlined  in  the  previous  section.  Once  the 

2 

smallest  real  root  Si*  of  Eq.  (61)  has  been  calculated,  the  modal 
displacement  functions  u^  and  w^  of  the  gravest  mode  can  be  found. 
These  modal  displacements  are 


and 


uQ  - (r/a)[l  - [(1  +v)Mq]/[(3  + v)0‘  - 2](r/a)2],  (64a) 


w0  - (2R/ a)  [( 1 - vK.2  - 2]“l 


X {l  - 2[<l  ♦ v)flJl/[(3  ♦ - 2](r/a):}. 


(64b) 


At  a nodal  circle,  the  normal  displacement  function  w^,  if  it  exists 
vanishes.  From  Eq.  (55),  it  is  seen  that  a nodal  circle  exists  at 


- a>J [i 3 + vK.‘  - :]/[:(!  + v)s;‘]^ 


(65) 


provided  that 


0 


2/(1  - v) 


(66) 


or  provided  that  the  lowest  resonant  angular  frequency 


is  such  that 


w0  < R'l|2E/[(l  - v)p]}\ 


(b7) 


Equations  (b6)  and  (67)  follow  from  Eq.  (65)  owing  to  the  fact  that 
the  quantity  (r/a)  cannot  exceed  unity. 
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SUMMARY 


In  this  report,  the  secular  equation  giving  the  resonant  frequencies 
of  a freely  vibrating  shallow  spherical  membrane  shell  has  been 
calculated  by  two  means.  The  exact  form  of  this  secular  equation,  under 
the  assumptions  of  the  theory,  is  given  by  Eq.  (AO).  An  approximate 
form  of  the  secular  equation  resulting  from  the  Ray leigh-Ritz  analysis 
is  given  by  Eq.  (61).  The  modal  displacement  functions  for  normal  and 
tangential  motions  that  are  associated  with  the  roots  of  the  exact 
secular  equation  are  respectively  given  by  Eqs.  (43)  and  (44). 
Corresponding  to  the  roots  of  the  approximate  secular  equation,  is  the 
expression  for  the  normal  modal  displacement  functions,  given  by  Eq. 
(64b),  and  the  expression  for  the  tangential  modal  displacement 
functions,  given  by  Eq.  (64a). 

REFERENCES 

[1]  H.  Kraus,  Thin  Elastic  Shells,  An  Introduction  to  the  Theoretical 
Foundations  and  the  Analysis  of  Their  Static  and  Dynamic  Behavior. 

(John  Wiley  & Sons,  Inc.,  New  York,  1967). 

[2]  C.  L.  Dym,  Introduction  to  the  Theory  of  Shells.  (Pergamon  Press, 

New  York,  1974)  pp.  8,  16. 

[3]  E.  Relssner,  "Stresses  and  Small  Displacements  of  Shallow 

Spherical  Shells",  J.  Math  Phys.,  2.5,  80-85,  279-300  (194b),  also 

27 . 240  (1948),  and  ^38,  lb-35  (1959).  These  references  were  obtained 
from  Ref.  [l],  p.  229. 


[4]  Ref.  [2],  pp.  21-22. 


